function fun_lorenz_Synchronous
%  Lorenz ????
%
% Author:LDYU
% Author's email: ustb03-07@yahoo.com.cn
%
% downloaded:
% http://ustbchaos.blogspot.com/2007/03/chaos-synchronization.html
 
[T,Y]=ode45(@lorenzf1,[0,3],[1;-10;1;10;2]);
subplot(2,1,1)
plot(T,Y(:,3),'k')
hold on
plot(T,Y(:,5),'b')
 
[T,Y]=ode45(@lorenzf2,[0,3],[10;10;10;-20;2]);
subplot(2,1,2)
plot(T,Y(:,3),'k')
hold on
plot(T,Y(:,5),'b')
 
%-------------------------------------
function dx=lorenzf1(t,x)
% dx=   a*(y-x)
% dy=  -x*z+r*x-y
% dz=   x*y-b*z
% dy1= -x*z1+r*x-y1
% dz1=  x*y1-b*z1
 
a=16;b=4;r=45.92;
dx=[ a*(x(2)-x(1))
    -x(1)*x(3)+r*x(1)-x(2)
     x(1)*x(2)-b*x(3)
    -x(1)*x(5)+r*x(1)-x(4)
     x(1)*x(4)-b*x(5)];
 
%-------------------------------------
function dx=lorenzf2(t,x)
% dx= a*(y-x)
% dy=-x*z+r*x-y
% dz= x*y-b*z
% dx1=a*(y-x1)
% dz1=x1*y-b*z1
 
a=16;b=4;r=45.92;
dx=[ a*(x(2)-x(1))
    -x(1)*x(3)+r*x(1)-x(2)
     x(1)*x(2)-b*x(3)
     a*(x(2)-x(4))
     x(4)*x(2)-b*x(5)];
 